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ABSTRACT 

The Cassini measurements of C 2 H 2 and C 2 H 6 at ~5 mbar provide a con¬ 
straint on meridional transport in the stratosphere of Jupiter. We performed a 
two-dimensional photochemical calculation coupled with mass transport due to 
vertical and meridional mixing. The modeled prohle of C 2 H 2 at latitudes less 
than 70° follows the latitude dependence of the solar insolation, while that of 
C 2 H 6 shows little latitude dependence, consistent with the measurements. In 
general, our model study suggests that the meridional transport timescale above 
5-10 mbar altitude level is >1000 years and the time could be as short as 10 
years below 10 mbar level, in order to £t the Cassini measurements. The derived 
meridional transport timescale above the 5 mbar level is a hundred times longer 
than that obtained from the spreading of gas-phase molecules deposited after the 
impact of Shoemaker-Levy 9 comet. There is no explanation at this time for this 
discrepancy. 

Subject headings: planetary systems—radiative transfer—atmospheric effects— 
planets and satellites: individnal (Jnpiter)— methods: nnmerical 


1. Introduction 

Meridional advection and mixing processes in the atmosphere of Jupiter are poorly 
known. Based on the Voyager infrared spectrometer data, several efforts to derive the at¬ 
mospheric circulation have been published (e.g., Gierasch et ah 1986; Conrath et ah 1990; 
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West et al. 1992). The first direct and quantitative derivation of meridional transport pro¬ 
cesses is based on the introduction of aerosol debris into the atmosphere of Jupiter by comet 
Shoemaker-Levy 9 (SL9). Friedson et ah (1999) conclude that the advection by the residual 
circulation calculated by West et al. (1992) is insufficient to explain the temporal distribution 
of cometary debris. Meridional eddy mixing coefficients on the order of 1-10x10^° cm^ s“^ 
are inferred in the regions between ~10 and 100 mbar. Later, based on the time evolution 
prohles of CO, CO 2 , CS, HCN, and H 2 O gas-phase molecules deposited after the SL9 impact, 
values of meridional eddy mixing coefficients as high as 2-5x10^^ cm^ s“^ are derived for 
pressures between ~0.1-0.5 mbar (Lellouch et al. 2002; Moreno et al. 2003; Griffith et al. 
2004). 

The Cassini measurements of stratospheric C 2 H 2 and C 2 He (Kunde et al. 2004) provide 
good tracers for characterizing mass transport in the upper atmosphere of Jupiter. How¬ 
ever, the weighting function of these observations is such that they are most sensitive to the 
altitude level near 5 mbar. Kunde et al. show that, at latitudes equatorward of ~70°, the 
relative magnitude of the abundance (or emission line intensity) of C 2 H 2 follows the latitu¬ 
dinal variation in solar insolation, while the abundance (or emission line intensity) of C 2 H 6 
is constant with latitude. Consequently, Kunde et al. conclude that the stratospheric merid¬ 
ional transport timescale at latitudes <70° derived from these Cassini data falls between the 
lifetimes of C 2 H 2 and C 2 H 6 . Due to the complexity in the auroral regions (contamination of 
line emissions from higher atmosphere due to temperature enhancement), we focus on the 
regions with latitude less than ~70° in this paper. 


2. Two-Dimensional Transport Model 

Because of the rapid rotation of Jupiter and its strong stratospheric zonal wind, the 
zonal variations in abundances should be minimized^. Therefore, a two-dimensional (2-D) 
model including net meridional and vertical transport should be an adequate hrst-order 
simulation of these Cassini observations of the trace hydrocarbon species. In the 2-D mode, 
the Caltech/JPL coupled chemistry/transport code, solves the mass continuity equation: 

dn^{y^z,t) ^ ^ = Pi{y, z,t) - Li{y, z,t), (1) 

where n* is the number density for the species i, ipi the transport flux, Pj the chemical 
production rate, and L* the chemical loss rate, all evaluated at time t, latitudinal distance 


^The rotation period and atmosphere-radiative timescale is ~10 hours and ~1000 days (at ~10 mbar) 
(Flasar 1989), respectively. 
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y and altitude 2 :. Reported herein are results for diurnally averaged steady-state solutions, 
i.e., <dni/dt> 0. 

For one-dimensional (1-D) problems {y dependence in eqn.(l) vanishes), the Caltech/JPL 
code (see, e.g., Gladstone et ah 1996) integrates the continuity equation including chemistry 
and vertical diffusion for each species using a matrix inversion method that allows large time 
steps. To take advantage of the computational efficiency of the 1-D solver, a “quasi 2-D” 
mode, which is a series of 1-D models at different latitudes coupled by meridional transport, 
has been developed. This quasi 2-D simulation has been tested against a case which has a 
known solution (Shia et al. 1990). Since the meridional transport in the stratosphere is not 
well understood, the quasi 2-D model uses a simple parameterization for mixing between 
the neighboring 1-D columns to simulate the meridional transport, i.e., Kyydrii/dy is added 
to the results of the 1-D computations, where Kyy is the meridional mixing coefficient. Our 
work provides an order of magnitude estimate of the meridional transport in the stratosphere. 

The vertical mixing coefficients {Kzz) and temperature prohle are taken from Glad¬ 
stone et al. (1996) and are assumed to be independent of latitude in the reference model. 
Sensitivities of the results due to the variations of temperature and K^z prohles are also 
presented. 

The meridional mixing Kyy will be determined by htting the Gassini measurements 
(Kunde et al. 2004). As a first order approximation, the Kyy prohle will be assumed also to 
be latitude-independent. 

The model atmosphere is gridded latitudinally and vertically in 10 and 131 layers, 
respectively. The vertical grid size is chosen to insure that there are >3 grid points in one 
pressure scale height (to achieve good numerical accuracy). The latitude grid points are at 
10, 30, 50, 70, and 85° in each hemisphere. 

The photochemical reactions {Pi and Li) are taken from Moses et al. (2000). At all 
latitudes, the mixing ratio of GH 4 in the deep atmosphere is prescribed to be 2.2x10“^, 
and the atomic hydrogen inhux from the top atmosphere is hxed at 4x10® cm“^ s“^ (see, 
e.g., Gladstone et al. 1996). In order to prevent the solution from oscillating seasonally, the 
inclination angle of Jupiter is prescribed to be zero (~3° actually). 


3. Simulation Results 

The Cassini measurements of the latitudinal distribution of C 2 H 2 and C 2 H 6 are repro¬ 
duced in Fig. 1; the contribution function for these is near 5 mbar (Kunde et al. 2004). 
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Since the absolute abundances of C 2 H 2 and C 2 H 6 have not been determined from these mea¬ 
surements, the results shown are arbitrarily scaled following the approach of Kunde et ah 
Because the C 2 H 2 abundances follow the latitudinal distribution of insolation (insolation is 
proportional to cosine of latitude), the meridional mixing timescale is expected to be longer 
than the C 2 H 2 photochemical lifetime. On the other hand, since the C 2 H 6 abundance ap¬ 
pears to be constant with latitude, the meridional mixing timescale is expected to be shorter 
than the C 2 H 6 photochemical lifetime. 

A more quantitative estimate for Kyy can be derived from the 2-D model. A validation 
of the 2-D model is from a 1-D model for Jovian hydrocarbon chemistry (Gladstone et al. 
1996; Moses et al. 2005), which reproduces extensive observations of hydrocarbon species 
as well as He 584 A and H Lyman-o; airglow emissions at low latitudes (observations are 
summarized in the Tables 1 and 5 of Gladstone et al. 1996 and Fig. 14 of Moses et al. 2005). 
The chemical loss timescales for G 2 H 2 and G 2 He drawn from our current work are shown 
in Fig. 2, along with vertical mixing timescales (also from our current work). At 5 mbar, 
the chemical loss timescales for G 2 H 2 and G 2 He are about 10 and 2000 years, respectively. 
Therefore, to reproduce the Gassini distributions of G 2 H 2 and G 2 H 6 , the meridional mixing 
time must fall within 10 and 2000 years at 5 mbar. 

Many simulations were performed with the quasi-2D model to explore the sensitivities 
of the abundances G 2 H 2 and G 2 H 6 at 5 mbar to the choice of Kyy. The results were param¬ 
eterized in terms of the ratio of the abundance at 70° latitude to the abundance at 10°. In 
general, the altitude variation of Kyy leading to a latitudinal gradient for G 2 H 2 consistent 
with the Gassini observations also led to a sharp reduction in the G 2 He abundance from 
equator to near-pole. Alternatively, for many cases tested, if the G 2 H 6 equator to near-pole 
variation was small, the same was true for the G 2 H 2 equator to near-pole variation. Shown 
in Table 1 are the model results that most adequately reproduce the Gassini observations. 
We found that a ‘transition’ level somewhere around 5 and 10 mbar must be present, in 
order to match the Gassini measurements. Above the transition altitude, Kyy <10® cm^ s“^. 
Below the transition altitude, Kyy > 2x10^® cm^ s“^, consistent with the analysis of the 
temporal spreading of the SL9 debris (Friedson et al. 1999). It is interesting to note that 
the results of Friedson et al. predict a reversal in the direction of the meridional component 
of the annual-mean residual velocity across the 5 mbar level in the southern hemisphere. 
This reversal represents a change with altitude in the relative contributions to the total 
annual-mean meridional heat flux from the component associated with the Eulerian-mean 
meridional velocity and the component associated with eddies. It is therefore possible that 
a change in the efficiency of meridional transport accompanies the reversal. 

These model results can be understood in the context of the basic photochemistry 



controlling the production and loss of C 2 H 2 and C 2 H 6 as outlined in Gladstone et al. (1996). 
In the atmosphere of Jupiter, most of hydrocarbon compounds are synthesized in the regions 
above ~0.1 mbar (shown in the shaded area in Fig. 2). Below 0.1 mbar, the production of 
C2H2 through CH4 photolysis mediated reactions is insufficient, because the process requires 
UV photons with wavelengths <130 nm, which has been self-shielded by CH4. In this region, 
the photolysis of C 2 He (<160 nm), which decreases strongly toward lower regions of the 
atmosphere, dominates the production of C 2 H 2 . The major chemical loss of C 2 H 2 and C 2 H 6 
is through hydrogenation and photolysis, respectively. As a consequence of this chemistry, 
C 2 H 2 is close to being in photochemical steady-state in the regions below ~5 mbar altitude 
level and its vertical gradient is large in this altitude range. Above 5 mbar level, transport 
is important for C 2 H 2 (Fig. 2). On the other hand, the abundance of C 2 H 6 is controlled 
by transport and its vertical gradient is very small (see, e.g.. Fig. 14 in Gladstone et al. 
1996). Therefore, if meridional mixing is sufficiently rapid below the transition altitude to 
uniformly mix G 2 H 6 with latitude, the tendency for uniform latitudinal mixing occurs also 
at 5 mbar. While above the transition altitude, uniform latitudinal mixing of G 2 H 6 also 
results in uniform latitudinal mixing of G 2 II 2 at 5 mbar. 

Fig. 3 shows the 2-D distribution of G 2 H 2 and G 2 lle calculated with the reference model, 
in which Kyy = 2x10^'^ below 5 mbar level and 2x10® cm^ s“^ above. Additional observations 
at different levels in the atmosphere can help constrain the 2-D dynamical properties of the 
Jovian atmosphere. 

vertical prohles of G 2 H 2 and G 2 H 6 at two latitudes. For comparison, 1-D results at low 
latitude are also overplotted. Additional observations at different levels in the atmosphere 
can help constrain the 2-D dynamical properties of the Jovian atmosphere. 

In the above calculations whose results are summarized in Table 1, assumptions were 
made with respect to temperature and vertical eddy diffusion coefficients independent of lat¬ 
itudes and in the selection of chemistry reaction parameters. Several additional calculations 
were performed to assess the sensitivity of the derived values for Kyy with respect to these 
assumptions. In one sensitivity test, the temperature prohle was progressively increased 
from 10° latitude so that, by 85°, the temperature prohle was 10% larger. Table 1, model 
G, shows that, with the adjusted temperature distribution, the Kyy that best simulated the 
latitudinal variations in G 2 II 2 and G 2 H 6 in the Gassini observations is the same as derived 
above. This is consistent with the conclusions of Moses and Greathouse (2005) who found 
little sensitivity to temperature in their calculations. In a similar fashion, Kzz was modihed 
linearly so that the value at 85° latitude was 10 times lower than the value at 10° latitude; 
enhanced Kzz at high latitudes cannot reproduce the measurements. The Kyy values that 
best simulated the Gassini observations (Table 1, model H) were again those derived earlier 
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in this paper. Finally, the same result for Kyy was found when the model chemistry was 
updated to be consistent with the reaction coefficients in Moses et ah (2005) (Table 1, model 
F). Therefore, the conclusions in this paper for the magnitude of meridional mixing as a func¬ 
tion of altitude are robust with regard to reasonable selection of atmospheric temperature, 
vertical mixing, and chemistry. 

4. Conclusion 

Our model simulation results have two implications. First, the meridional transport 
time as short as 10 years {Kyy ^ 10 ^^ cm^ s“^) exists only in the altitude range below the 
10 mbar level. Second, above a transition level somewhere between 5 and 10 mbar, the 
meridional transport time is not shorter than 1000 years {Kyy < 10® cm^ s ^). While these 
inferred Kyy values for the atmosphere at and below the 5-10 mbar level are consistent with 
the conclusions of Friedson et al. (1999) derived from an analysis of the SL9 debris evolution 
with time, the Kyy value above the transition level is much smaller than that ( 10 ^^ cm^ 
s“^) derived from analysis of the time evolution of the distributions of gas-phase trace species 
deposited after the SL9 impact (Lellouch et al. 2002; Moreno et al. 2003; Griffith et al. 2004). 
There is no explanation at this time for this discrepancy. 

It has been shown that CH 4 and C 2 He contribute to the heating and cooling of the 
stratosphere of Jupiter, respectively (Yelle et al. 2001). To a hrst order approximation, 
the cooling at/below 5 mbar level would be constant with latitude, being determined by 
two factors, the temperature and the abundance of C 2 He. These two are nearly constant 
between the equator and mid-latitudes (Flasar et al. 2004; Kunde et al. 2004, also Fig. 3). 
The heating function, however, is sensitive to the magnitude of the solar insolation and, 
consequently, it is a good approximation to assume that this function has the same latitude 
dependence as the solar insolation. Therefore, the global circulation could be driven by the 
differential heating between latitudes. This circulation driven by heating through absorption 
of radiation by gas-phase molecules will provide a hrst order estimate of the importance of 
aerosol heating in the stratosphere. 

This research was supported in part by NASA grant NAG5-6263 to the Galifornia In¬ 
stitute of Technology. Special thank to Julie Moses for her insightful comments. 
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Fig. 1.— Relative abundances of C 2 H 2 (upper) and C 2 H 6 (lower) at 5 mbar as a function 
of latitude. Asterisks are the Cassini measurements (Kunde et al. 2004). Diamonds and 
circles are, respectively, high-latitude Cassini measurements that do include (diamonds) and 
do not include (circles) auroral longitudes. The calculated abundances of C 2 H 2 and C 2 H 6 
are normalized to those at the equator. Solid line represents model results with the reference 
Kyy\ constant 2x10^'^ cm^ s“^ below the 5 mbar altitude level and 2x10® cm^ s“^ above 
(model B). Dotted line represents the cosine function of latitude. 
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Fig. 2.— Timescales for the chemical loss of C 2 H 2 (solid lines) and C 2 H 6 (dashed lines) and 
for vertical transport (dotted line). The vertical transport timescale is dehned by 
where K^z and H are the vertical diffusion coefficients of CH4 and atmospheric scale height, 
respectively; values of time constants are derived from our reference model. Thick and thin 
lines represent values at latitudes 10° and 70°, respectively. The horizontal arrow indicates 
5 mbar level where peak in the contribution function for the Cassini measurements of C 2 H 2 
and C 2 H 6 lies. The vertical long-dashed line is a meridional mixing time equal to R^j/Kyy, 
where Rj is the radius of Jupiter and Kyy = 2x10^*^ cm^ s“^. The shaded area shows the 
photochemical production region for the hydrocarbons. 







Table 1. Summary of model results 



Chemistry^ 

Temperature^ 

K ^ 

Transition*^ 

Kyy below° 

Kyy above° 

C2H2‘1 

C2H6'^ 

Solar flux 







0.34 

0.34 

Cassini 







0.50 

0.87 

Model A 

Standard 

Standard 

Standard 

5 

2x10^° 

0 

0.49 

0.77 

Model B 

Standard 

Standard 

Standard 

5 

2xl0l0 

2x10® 

0.54 

0.80 

Model C 

Standard 

Standard 

Standard 

10 

2x10^° 

2x10® 

0.46 

0.73 

Model D 

Standard 

Standard 

Standard 

10 

2xl0ii 

0 

0.46 

0.75 

Model E 

Standard 

Standard 

Standard 

10 

2xl0ii 

2x10® 

0.51 

0.79 

Model F 

Moses et al. (2005) 

Standard 

Standard 

5 

2xlOio 

2x10® 

0.59 

0.84 

Model G 

Standard 

xl.l 

Standard 

5 

2x10^° 

2x10® 

0.49 

0.76 

Model H 

Standard 

Standard 

xO.l 

10 

2xl0ii 

2x10® 

0.48 

0.91 


^ Standard chemistry is taken from Moses et al. (2000). Standard temperature and Kzz profiles are from Gladstone et al. 
(1996). Modified chemistry is from Moses et al. (2005). Reference temperature and Kzz profiles are set at equator (10°). 
The maximum changes in temperature (increase by 10%) and Kzz (reduced by a factor of 10) profiles are at 85°, with the 
change assumed to be linearly proportional to the angle of latitudes. 

^ Transition level of Kyy profile. The level is given in units of mbar. 

° Values of Kyy below and above the transition level. Kyy is in units of cm^ s~^. 

^ Ratio of abundance at ±70° latitude to abundance at the equator. Values for Cassini measurements are averaged at 
±70°. 







